(12) INTERNAffONAL APPUCATION PUBLISHED UNDER THE PATENT COOPERATION TREATY (PCT) 



(19) World Intellectual Property OrganizatioD 
Intemational Bureau 




(43) International Publication Date (10) Intemational Publication Number 

30 November 2000 (30.11.2000) pCT WO 00/71980 Al 



(51) International Patent ClassificatioD^: G01J3/28, 

C301N 21/27 

(21) Intemational AppUcation Number: PCT/CA00/005S3 

(22) Intemational Filing Date: 23 May 2000 (23.05.2000) 

(25) filing LangD^e: English 

(26) Publication Language: English 

(30) Priority Data: 

2^73,132 21 May 1999 (21.05.1999) CA 

(71) Applicant (for all designated States except US): MEA- 
SUREMENT MICROSYSTEMS A-Z INC. [CA/CA]; 
1 14 Place des Chenaux, Trois-Rivi^, Qjj^aec G8 Y 5W5 
(CA). 



(72) Inventors; and 

(75) Inventors/AppUcants (ft>r US only): BARVHCZ, An- 

drzej [CA/CA]; 114 Place des Chenaux, Trois-Rivi^res, 
Quebec G8Y 5W5 (CA). \VISN1£VVSKT, Michael,. P. 
[PL/CA]; 1125, boul. Des Recollets, App. 23. Trois-Riv- 
ieres, (Quebec OaZ 3X3 (CA). MORAWSKI, Roman, 
Z. IPL/PL]; Filtrowa 71 M8. PL^2-055 Warsaw (PL). 
SLIMA, Mohamed, Ben [TN/CA]; 4025 Louis Pinard 
ApL #2, Trois-Rivi^, Quebec G8Y 4C2 (CA). 

(74) Agent: FREEDMAN, Gordon; Fieedman Sc Associates, 
117 Centrepolnte Drive, Suite 350, Nepean, Ontario K2G 
5X3 (CA). 

(81) Designated States (nathnal): AE, AG, AL, AM, AT, AU, 
AZ, BA, BB, BG, BR, BY, CA, CH. CN, CR, CU, CZ, DE, 
DK, DM, DZ, EE, ES, H, GB, GD, GE, GH, GM, HR, HU, 
ID, IL, IN, IS, JP, KE, KG, KP, KR, KZ, LC, LK, LR, LS, 
LT, LU, LV, MA. MD. MG. MK, MN. MW, MX, MZ, NO, 

[Continued on next page] 



(54) Tide: METHOD AND SYSTEM FOR ADAPTIVE INTERPRETATION OF SPECTROMETRIC DATA COMBINED WITH 
CONTINUAL RE-CAUBRAnON 



CAUBRATION OP THE 
SPECTROMETRIC DEVICE 
(sub pfoeedore ARICRC cal) 




Spectra of 
■ caUbrmtessBpIci 



Spcdnof 
calibntiossam 



re<:aubration of the 
spectrometric device 

(sob jiroccdore ARICRC.ncal) 



ADAPTIVE RECONSTRUCTION 

OF HIE SPECTRUM 
(Mb jncidun ARlCRC^ree) 



> 



I 



ADAPTIVE INTERPRETATION 

OFTUESPECraUM 
(sQb_proccd«r« ARXCRC.fnt) 



(57) Abstract: A method 
for adaptive reconstniction 
and inteipretation of spectra 
combined witb ie*caIibration of 
a device providing spectrometric 
data, consist in canying out an 
automatic calibration using the 
external reference ligbt spectrum 
and its corresponding digital 
reference data stored in the internal 
memory of the device. The 
continual le-calibration procedure 
allows for automatic ads^tion 
of the values of the coefficients 
in a reconstruction sub-procedure 
as well as the estimation of the 
values of the coeffideats in an 
interpr^ation sub-procedure on 
the bases of a current shape of 
spectrometric data. 



wo 00/71980 Al liHi 



niiiiiiHiiiiiiiin 



NZ. PL, FT, RO. RU, SD, SE. SG, SI. SK, SL, TJ, TM, TR, 
TT. TZ. UA. UG, US, UZ. VN, YU. ZA, ZW. 

(84) Designated States (regionaiy. ARIPO patent (GH. GM. 
KE, LS, MW. MZ. SD. SL. SZ. TZ, UG, ZW), Eurasian 
patent (AM, AZ, BY, KG. KZ. MD. RU. TJ, TM). European 
patent (AT. BE, CH. CY. DE, DK, ES , H. FR, GB . GR, IE, 
rr, LU. MC. ML. PT. SE). OAPI patent (BP. BJ, CP, CG, 
a. CM, GA. GN. GW. ML. MR. NE, SN, TD, TG). 



Published: 

— With international search report. 

— Before the expiration of the time limit for amending the 
claims and to he republished in the event cf receipt of 
amendments. 

For two-letter codes and other abbreviations, refer to the "Guid- 
ance Notes on Codes and Abbreviations" appearing at the begin- 
ning of each regular issue of the PCT Gazette. 



wo 00/71980 



PCT/CA00/00S83 



Method and system for adaptive interpretation of spectrometric data 
combined with continual re-calibration 

Field of Invention 

This invention relates generally to light-spectrum-measurement and niorc 
5 specifically to compact systems for light-spectrum-measurement for portability and for //; 
situ applications. 

Background of the invention 

Spectroscopy is an analytic technique concerned with the measurement 
characterization of the interaction of radiant energy with matter, including the instruments 

10 designed for this purpose, called spectrometers ^ and corresponding means of the 

interpretation of the interaction both at the fundamental level and for practical analysis. 
The distribution of radiant energy, absorbed or emitted by a sample of a substance under 
study, is called its spectrum. Interpretation of spectra provides fundamental information 
at atomic and molecular energy levels, the distribution of species within those levels, the 

15 nature of processes involving change from one level to another, molecular geometries, 
chemical bonding, and interaction of molecules in solution. At the practical level, 
comparisons of spectra provide a basis for the determination of qualitative chemical 
composition and chemical structure, and for quantitative chemical analysis. 

The spectra to be measured have two conmion features: they are non-negative, 
20 and they may be decomposed into relatively flat regions and peaks of various widths. The 
interpretation of the data representative of those spectra obtained by means of 
spectrometers, very frequently involves estimation of some parameters of the peaks, viz.: 
of their positions and magnitudes (i.e. areas or heights). The parameters of the peaks are 
used, in particular, for qualitative and quantitative analysis of complex chemical 
25 substances: 

• the positions of peaks are used for identification of the compounds which are 
present in a sample of the substance under study; 
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• the magnitudes of peaks are used for estimation of the concentrations of the 
identified compoimds. 

The estimation of the parameters of spectrometric peaks directly from the 
acquired spectrometric data is in many practically important cases problematic because 
5 of: 

• instrumental imperfections of the spectrometer, and 

• some natural phenomena in a sample of the substance under study. 

The imperfections of the spectrometer produce blurring of peaks and noise-type 
disturbances in the data representative of the measured spectrum. On the other hand, 

1 0 some quantum phenomena in the substance enlarge considerably the width of peaks. Both 
instrumental blurring and natural enlargement of peaks may contribute to their 
overlapping and fusion, making impossible their proper identification. Different 
algorithms of deconvolution have been developed for reducing instrumental blurring; a 
review of these algorithms may be found in Jansson P. A. (Ed.) : Deconvolution with 

15 Applications in Spectroscopy. Academic Press, 1984 and Morawski R. Z., 

Szczecin ski L., Barwicz A.: "Deconvolution Algorithms for Instrumental Applications - 
A Comparative Study", J. Chemometrics, 1995, VoL 9, pp. 3-20. These algorithms are 
designed and assessed using the criteria of the quality of deconvolution, which is not 
specific of spectrometry; e.g. the root-mean-square errors of approximation. 

20 Consequently, the estimates of the magnitudes determined on the basis of the results of 
deconvolution may turn out to be very poor. On the other hand, the curve-fitting 
algorithms may be very inefficient if not provided with a good initial guess of the sought- 
for estimates. These two observations inspired independently many researchers to 
subsequently use a deconvolution algorithm for estimation of the positions and a curve- 

25 fitting algorithm for estimation of magnitudes of spectral peaks. 

Both the algorithms of spectrum reconstruction and the algorithms of spectrum 
parameter estimation require considerable amount of information on the mathematical 
model of the spectrometric data. This information is acquired during initial calibration of 
the spectrometer. 
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The numerical algorithms necessary for calibration of the spectrometer, 
reconstruction of the spectrum and estimation of its parameters are available in many 
software libraries, e.g. in some MATLAB Toolboxes, including specialized libraries 
dedicated to processing of spectrometric data, e.g. GRAMS/32 of Galactic Industries. In 
5 the patent domain, there is only a partial methodology for spectrum deconvolution, e.g. 
the US Patent No.5 247 175, issued September 21, 1993, in the name of 
Alan E. Schoen et aL^ entitled Method and Apparatus for the Deconvolution of 
Unresolved Data as well the US Patent No.4 941 101, issued July 10, 1990, in the name 
of Paul B. Crilly, entitled Method for Analyzing Chromatograms specific to the domain of 
10 chromatography. On the other hand, there is a good example of a complete methodology 
of spectrometric data interpretation - c.f the US Patent No. 5,991,023, issued November 
23, 1999, in the name of R. Z Morawski et al., entitled Method of interpreting 
spectrometric data. 

However, in all the above-listed references, none of them offers full support for 

15 the methodology of on-line continual re-calibration and adaptive spectrometric data 
interpretation, mainly due to the fact that there are no closed procedures of such re- 
calibration and adaptive interpretation. 

The demand for an increased-resolution and a self calibrated spectrometric device 
comes from everyday practice. The US Patent No. 5040889, issued August 2 1 , 1 99 1 , in 

20 the name of Thomas J. Keane, entitled Spectrometer with Combined Visible and 

Ultraviolet Sample Illumination^ describes a spectrometer with an automatic calibration 
of tiie photodetectors, based on the use of a pivotable standard white sample positioned in 
the path of the light and controlled by a computer system. A calibration factor is 
computed and recalculated for each photodetector whenever the temperature changes 

25 within the spectrometer, and then used to multiply the output signal of the photodetectors. 
However, this is a partial static calibration, since it corrects only the effect of the 
photodetectors and doesn't take into account the effect of the dynamical errors introduced 
by the light diffraction element and the slits. Thus, there is a need for another dynamic 
calibration procedure to eliminate these instrumental errors and enhance the spectrum 

30 resolution. In the fields of telecommimications, especially for optical performance 
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monitoring of networks, there is a need for a spectrometric device, which must be 
compact, operating under a broader range of environmental conditions and more 
particularly self-calibrating. In Miller C, Pelz L.: "Fabry-Perot Tunable Filters Improve 
Optical Channel Analyzer Performance", Lightwave, March 1999, Vol. 16, No. 3, a card 

5 mounted optical channel analyzer performance, based on the use of a fiber Fabry-Perot 
tunable filter, is described. This device uses an optical switch to alternately scan the light 
spectrum (coming for example from a DWDM system) and a known wavelengths 
reference source. This latter is used for the continual calibration of the device in order to 
correct the thermal drifts, the nonlinearities of the actuator and the voltage variations. A 

10 deconvolution algorithm is also used to remove the effects of the tunable filter spectral 
response in order to enhance the resolution of the measured light spectnim. The major 
drawback of this device is that it measures one wavelength or path length difference at a 
time, thus necessitating a long time to process all flie optical channels. Thus, a need 
remains for a low-cost miniaturized spectrometric device, self-calibrated, with a spectral 

15 resolution comparable to that of conventional optical spectrum analyzer, and capable of 
determining the spectral signature of a wide variety of light spectra in sifu. 

In the domain of spectroscopy, the demand for an adaptive recoastruction and 
interpretation procedures combined with a continual re-calibration is generated by the 
need for a spectrometric device that can be compact, autonomous with a higher degree of 
20 reliability, self-calibrated, operating under a broader range of environmental conditions, 
and be cost-effective for wide deployment. The main objective of the invention is to 
respond to this demand. 

Object of Invention 

It is an object of this invention to provide a method for adaptive interpretation of 
25 spectra on the basis of spectrometric data. 

It is an object of this invention to provide a method for calibration and continual 
re-calibration of integrated spectrometric devices allowing for miniaturization and low- 
cost mass production of micro-optical spectrum analyzers. 
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Summary of the invention 

It is an object of this invention to provide a procedure for interpretation of 
spectrometric data, representative of a sample of the substance under study, comprising 
the steps of: 

5 • initial calibration of the spectrometric device; 

• continual re-calibration of the spectrometric device 

• adaptive reconstruction of the spectrum under study; and 

• adaptive interpretation of the spectrum under study. 

The invention also provides a method for calibration and continual re-calibration 
10 of the spectrometric device, i.e. for identification of a numerical relation between the 
measurement results produced by this spectrometric device and established measurement 
standards of spectrum. 

The invention also provides a method for continual re-calibration of the 
spectrometric device allowing for automatic adaptation of the method to possible 
15 variation of the instrumental imperfections and a shape of spectrometric data due to the 
aging of the device and/or thermal drifts and nonlinearities of the detector. 

The invention also provides a method for numerical adaptive reconstruction of the 
spectrum under study, combined with the continual re-calibration of the spectrometric 
device, i.e. for partial elimination of the effects of noise and blur corrupting the 
20 spectrometric data. 

The invention also provides a method for numerical adaptive interpretation, i.e. 
estimation of the positions and magnitudes of peaks the spectrum of the light is composed 
of. 

According to the invention, the calibration of the spectrometer is accomplished 
25 using an auxiliary light whose spectrum is assxuned to be known and representative of the 
spectra to be measured in the series of experiments following the calibration. The result 
of calibration has the form of two operators defining the following numerical algorithms: 
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• an algorithm for simulation of the spectrometric data, given the spectrum to 
be measured (operator of projection 

• an algorithm for estimation of the spectrum, given the data (operator of 
reconstruction 91). 

5 The first algorithm is an algorithm of model identification including choice or 

estimation of the structural parameters of the model and estimation of its functional 
parameters. The second algorithm is an algorithm of deconvolution, or an algorithm of 
generalized deconvolution, or an algorithm for numerically solving the first-kind integral 
equations. 

10 According to the invention, the continual re-calibration of the spectrometer is 

accomplished using an external reference light spectrum and the digital reference data 
representing this reference light, stored in a computer or an internal memory, when 
available, of the spectrometric device. The re-calibration of the device is aimed at 
updating the parameters of the above mentioned operators. 

15 According to the invention, the processing of the data representative of the light 

under study is comprising two main operations: 

• estimation of the unknown spectrum or its idealized version on the basis of 
the acquired spectrometric data, using the operator of reconstruction S?; 

• estimation of the parameters of this spectrum on the basis of the acquired 
20 spectrometric data, using the operator of projection ^ and an algorithm of curve 

fitting. 

The method of adaptive interpretation of spectrometric data is usefiil for providing the 
quality light-spectrum measurement especially in the case of miniature integrated 
intelligent spectrometric devices. When combined with automatic recalibration 
25 continually or at intervals the spectrometric device keeps optimum performance, in 
particular, for teleconmnmication applications. 

The proposed method allows for automatic calibration of especially miniature 
integrated spectrometric devices during a manufacturing process and absent physical 
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tuning of individual devices. The calibration will compensate for many forais of 
fabrication discrepancies without modification of the physical devices. 

As tiie consequence, it contributes to a low-cost flexible quality light-spectrum- 
measurement solution, especially on-chip light spectrum measurement, which may be 
5 easily adapted to very diversified applications via customization of the specific set of 
algorithms. 

Brief description of the drawings 

Exemplary embodiment of the invention will now be described in conjunction 
with the drawings in which: 
1 0 Fig. 1 • presents an exemplary measuring system; 

Fig 2. presents the flow diagram of the procedures of adaptive reconstruction and 
interpretation 

a) the procedure AI-RC, 

b) the sub-procedure AI-RC_cal, 
15 c) the sub-procedure AI-RC_recal, 

d) the sub-procedure AI-RC_rec, 

e) the sub-procedure AI-RC_int; 

Fig. 3. presents the exemplary measuring system; 
Fig. 4. presents the spectrum of the test sample: 
20 a) real spectrum x{X) , 

b) the test data {y„}; 

Fig. 5. presents the spectrum of the sample used for calibration: 

a) real spectrum a), 

b) the data {^r'); 

25 Fig. 6. presents the output of: 

a) the rational filter, 

b) the spline-based Kalman filter; 

Fig. 7. presents the final result of spectrometric data interpretation; 
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Fig. 8. presents the exemplary measuring system - MM microOSA; 
Fig. 9. shows possible applications of MM microOSA™ in Dense Wavelength Division 
Multiplexing (DWDM) system; 
Fig. 10. presents the spectrum of the test light: 
5 a) the idealized spectrum s{X) of the light at the input of the DWDM transmitter, 

b) the test data {y„)\ 

Fig. 1 1. presents the emission spectrum of the laser; and, 

Fig. 12. illustrates the final result of adaptive reconstruction and interpretation of 

telecommunication data. 

10 Detailed Description of the Invention 

The procedure of adaptive reconstruction combined with continual re-calibration 
of a device providing spectrometric data and adaptive interpretation of spectra (AI-RC) is 
designed for a measuring system, for example the system shown in Fig.l, comprising: 

• a spectrometric device in the form of a spectrometer or its essential part, in the form 

15 of a micro-optical spectrum analyzer converting an optical signal carrying the information 
on the measured spectrum into a digital code representing this spectrum; 

• a processing means for processing the digital representation, in the form of a general- 
puq50se computer, a microprocessor, a general-purpose digital signal processor, or an • 
application-specific digital signal processor; and 

20 • other functional elements necessary for measuring the spectrum of light 

The following notation is used for the description of the procedure AI-RC: 
X - wavelength; A e[A^^,\,tix]* 
N - number of data acquired by a spectrometric apparatus; 
M - step of wavelength discretization; AA = - )/{N - 1) ; 
25 X„ - n-th datum acquired by the spectrometric apparatus; = A^^j^ +{n- l)AA for 
w = l,...,i^; 

x(A) - real spectrum of a light under study; 
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1 - a vector of positions of peaks within the spectrum x(A,) is composed of, 
i - an estimate of 1 ; 

a - a vector of magnitudes of peaks within the spectrum x(X) is composed of; 

5 a=[a, ^2 ... a^Y; 

a - an estimate of a; 

5(A;l,a) - an idealized spectrum of light under study and, in the preferred embodiment, 
assumed to have the form: 

K 

^(A;l,a) = X^jtV,(^/jt) 
10 where v,(A,/) is an isolated, normalized peak in 5(A;l,a), whose maximum is located at 

|^v,(A,/)JA = l for /€[^i„,A,3,]; 

{yn} " spectrometric data representative of x(X) , acquired by means of the spectrometric 
apparatus; 

15 [y„}^{y„\n = \,...,N]', 

x"'{X) - a real spectrum of light used for calibration of the spectroinetric device; 
s(A'X°^,z^^^ ) - an idealized spectrum of light used for calibration of the spectrometric 
device; 

} - the spectrometric data, representative of x'^(X) used for calibration of the 
20 spectrometric device; 

^recai ^ ^^^^ spcctTum of an external reference light used for the re-calibration 
of the spectrometric device; 
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}^reeal | ^ spectrometric data measured using the spectrometric device, 

representative of x^^^\X) used for the re-calibration of this device; 

^(A;r^^^^a'*^^^^) - an idealized spectrum of a light used for the re-calibration of the 

5 spectrometric device; 

G - an operator (transform) of projection mapping the idealized spectrum s{X\\^2i) mto 

the space of the data: 

{xj=e [5(A;l,a);p, ] 

where p^ is a vector or matrix of the parameters of the operator G , to be determined 
10 during calibration of the spectrometric device; p^ -\pq ,2 ^r: 



Pi5 = 



PG ,1,1 P& X2 
PGX\ PGX2 



R - an operator (transform) of reconstruction and interpretation: 
R =R 'oR ' 

R ' - an operator (transform) of reconstruction, in the form of a generalized 
15 deconvolution method, for transforming the data into an esiimate {x„} of x{a) : 

where Pr r = [p,^ r , p,^ r j are parameters of the operator R ^ including regularization 

parameters for determination during calibration of the spectrometric device. 
R * - an operator (transform) of interpretation, in the form of a generalized deconvolution 
20 method, for transforming the data {x„ } into an estimate S{X) of s(A; 1, a) : 

where 1 = 1 , 1 2 ^® parameters of the operator R * including regularization 
parameters for determination during calibration. 

2. Functional specification of a preferred embodiment of the invention 
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2.1. General specification of the procedure Al-RC 

The main objective of the procedure AI-RC is estimation of the positions i and 
magnitudes a of the peaks contained in a spectrum of a light under study x(X) on the 

basis of the acquired spectrometric data {y„] . The feasibility of this operation is 

5 critically conditioned by an auxiliary operation on the reference data {y } and 

corresponding reference spectrum x^°\X) , referred to as calibration of the spectrometric 
device. This operation is aimed at the acquisition of information relating to a 
mathematical model of a relationship between spectrometric data and idealized 
spectrometric data known to be measured, which underlies the method chosen for 
10 estimation of the parameters 1 and a. Although calibration does not necessarily 

immediately precede processing of each sequence or sample of spectrometric data (j/^j , 
preferably valid calibration results are always available during spectral processing. 

The main difficulty relating to estimation of the positions 1 and magnitudes a of 
spectrometric peaks within a sensed spectrum is implied by blurring of those peaks, 
15 caused by physical phenomena in a sample, and by blurring of their representations in the 

data [y^] , caused by an imperfect spectrometric apparatus. This difficulty is overcome in 

the procedure AI-RC through application of a series of steps for adaptive reconstruction 
and adaptive interpretation of an idealized - hypothetical - spectrum s(X'Xz) in order to 
correct the sensed data to reduce both types of blurring; if s{A; 1, a) is assumed to be an 
20 approximation of x{X) , then only the instrumental blurring needs to be corrected. 

In accordance with the above general fimctional specification, AI-RC is composed 
of the following broadly defined steps: 

• initial calibration of the device providing the spectrometric data, i.e. of the 
spectrometric apparatus or of the spectrometric sensor or of the spectrometric transducer 

25 (the sub-procedure AI-RC_cal), 

• continual re-calibration of that device on the basis of auxiliary reference data obtained 
by means of the device stimulated with the reference light (the sub-procedure AI- 
RC_recal), 

SUBSTITUTE SHEET (R ULE 26) 
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• adaptive reconstruction of the spectrum under study {x„}, i.e. estimation of this 
spectrum on the basis of the spectrometric data {j;} representative of the spectrum {xj 
(the sub-procedure AI-RC_rec), and 

• adaptive interpretation of the spectrum under study {x„ } in the form of adaptive 

5 estimation of the positions /, , , ... and magnitudes aj , , ... of peaks this spectrum is 
composed of - on the basis of an estimate {x„} of the real spectrum, obtained by means of 
the sub-procedure AI-RC_rec (the sub-procedure AI-RC_int). 

2.2. Detailed specification of the sub-procedure AI-RC cal 

10 The sub-procedure AI-RC_cal comprises the following operations: 

a) choosing the method for idealization of the spectrum, Le. the form of peaks of the 
idealized spectrum {s„} ; 

b) choosing the form of the operator R for reconstruction of the idealized spectrum: 
{f^} =R [{y„}; ], where p,^ is a vector of the parameters; 

15 c) choosing the rules for adaptation of the parameters p^^ to the selected features 
Fi, Fj, ... of the spectrum under study {x„ } or of the data {yj - on the basis of the 
library of parameters | y =1,...,/} to be determined during initial calibration, and 
updated during continual re-calibration; 

d) choosing reference light signals for the initial calibration: (x^'^O*)} and related 
20 {yfO0}fory=l,....J; 

acquiring the data ^^"^l representative of (x^'^'O)} J-U...,J ; 

e) pre-processing of the data ^^''} including the elimination of outliers, baseline 
subtraction, smoothing, acquiring a priori information for the operation 2.g (e,g, a pre- 
estimate of the variance of errors in the calibration data), normalization, e/c; 

25 f) determination of the vectors of parameters p^^ - on the basis of the triplets: 

{xfOOjand trO)}fory=l,...,J; 
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g) approximation of the relationship between the parameters P|^ and the features 
Fi , » *•* controlling the adaptation. 

2.3. Detailed specification of the sub-procedure AI«RC recal 

The sub-procedure AI-RC_recal comprises the following operations: 
5 a) choosing the reference light signals for re-calibraticn: (vr^'^'fy)} 2nd related {y^^'^'O')} 
for y=l,,..,J; 

b) acquiring the data ^7''''' j representative of {x^'^^'^C/)} for 7=1,..., J ; 

c) pre-processing of the data {y^^^^'j including elimination of outliers, baseline 

subtraction, smoothing, acquiring a priori information for the operation 3.d {e.g. a pre- 
10 estimate of the variance of errors in the calibration data), normalization, etc.\ 

d) updating of the vectors of parameters p^^ on the basis of triplets: ^T^iff^ 

and{yr'''t/")}for;=l,....J; 

e) updating the approximation of the relationship bet^^veen the parameters pp^ and the 
features ^1,^2,... controlling the adaptation. 

15 2.4. Detailed specification of the sub-procedure AT-RC rec 

The sub-procedure AI-RC_rec comprises the following operations: 

a) acquiring data |>?„| representative of the spectrum under study {x^]\ 

b) pre-processing of } including elimination of outiiers, baseline subtraction, 

smoothing, acquiring a priori information for the estimation of peak parameters in the 
20 form of a pre-estimate of the variance of errors in the data, acquiring a priori information 
for the reconstruction in the form of determination of apparatus fabrication imperfections 

on the basis of characteristics of measured spectrum |5^„ | , ambient temperature sensed 

with a temperature sensor, and aging of working apparatus elements, normalization, etc.; 

c) estimation of the spectrum under study {x^ }on the basis of the data {>;„} , by means of 
25 the operator R ' with the parameters , being continuously adapted to the features of 
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3. An exemplary embodiment of the invention 
Example 1 

A particular version of the procedure AI-RC has been designed for a measuring 
system shown in Fig. 3 and including: the absorption mini spectrophotometer - model 
5 S 1 000 by Ocean Optics®; and a personal computer PC. 

The following measurement parameters were selected both for calibration and for 
acquisition of the test data: 

• the wavelength range: = 450/im, = 675nm ; 

• the number of data acquired by the spectrophotometer: N = 1 024 ; and, 

10 •the step of wavelength discretization: = (A^ - >l^„ )/{N - 1) = 0.22nm . 

The test data were acquired for a standard hohnium perchlorate solution sample; 
its real spectrum x{A) is shown in Fig. 4a. The parameters of this spectrum are the 
following: 

• the vector of the positions of peaks: 

15 

I = [452.2 468.2 473.1 485 537.3 541.3 543.7 639.8 644.2 652.2 656.4f 

• the vector of the magnitudes of peaks: 

a=[l.l7 0.19 0.2 0.461 0.988 0.41 0.297 0.821 0.406 0.246 0.263f 
The idealized spectrum of a sample under study has been assumed to have the form: 

20 ^(A;I,a) = 2;^*v,(A,/,) 



with the peaks defined by: 



v,aJ)^S{X-l) for/€[4,j„,^3xl 



The set of data representative of x{A) , acquired by means of the spectrophotometer, 




SUBSTITUTE SHEET (RULE 26) 



wo 00/71980 PCT/CA00/00S83 

15 



Calibration data have been acquired for the standard holmium oxide filter sample; 
its real spectrum x^'^^iX) is shown in Fig. 5a. The parameters of this spectrum are the 
following: 

• the vector of the positions of peaks: 

5 

r'= [450.7 454.5 460.4 463.5 473.9 483.9 4S8.S 536.4 547.5 633.4 636.5 64S.4f 



9 



the vector of the magnitudes of peaks: 



10 a'"' =[0.486 0.799 0.949 0.58 0.125 0.117 0.114 0.284 0.106 0.129 0.149 0.1097 

. The idealized spectrum of a sample used for calibration si^iV"^ ^sl"""') has been assumed 

to have the form: 

.(A;l«-,a"')=i<"'v,(^/r') 

kml 

The set of data representative of x*^(X) , acquired by means of the spectrophotometer, 
15 {y„"'}s I n = 1,...,1024}, is shown in Fig. 5b. 

The chosen operator (transform) of projection, mapping the idealized spectrum 
i(A;I,a) into the space of the data: 

{>>„}=G[^(A;l,a);ps] 

has been defined by the following operations: 

20 x{X) = ex^j^gJX-A')ln[s{A'-Xn)]dX'^ 

y„'=){^,) for n = l,...,N 
The function g^yiX) has been assiimed to have the form of the Gauss function: 
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Consequently, the vector of the parameters of the operator G contains the discrete 
values of gj^d^) and the parameter o^. 

Tlie chosen operator (algorithm) of reconstrucLion, transforming the data [y^j 
into an estimate s{X) of 5(/l;l,a), 

5 iW=f<[{y.};Pft] 

is defined by in the following way: 

a discrete estimate |jc„| of is found by means of the rational filter described 

by SzczeciA ski L.: "M6thodes non lineaires de reconstitution de signaux pour des 

applications spectrometriques", Ph. D. thesis, Universite du Quebec - Institut National de 
10 la Recherche Scientifique «Teleconunimications», Montreal 1997. and applied to the data 

(?.}: 

• an estimate of 5(A;l,a) is computed by means of the adaptive rational filter as 
described in Wisniewski M. P., Morawski P^. Z., Eanvicz A., 1999, "An Adaptive 
Rational Filter For Interpretation of Spectrometric Data", to be published in IEEE Trmis, 
15 Instrum. & Meas. applied to (-^fl) - 

The vector Pr r = [p^ r | r ^ — f ^® parameters of the operator R contains 
the coefficients of the rational filter, as well as the discrete values of the fimction ^^(A) 
as described in Ben Slima M., Szczecinski L., Massicotte D., Morawski R. Z., 
Barwicz A.: "Algorithmic Specification of a Specialized Processor for Spectrometric 
20 Applictions", Proc. lEEEInstrum. & Meas. Technology Conf. (Ottawa, Canada, May 19- 
21, 1997), pp. 90-95 and in Ben Slima M., Morawski R. Z*, Barwicz A,: "Kahnan-filter- 
based Algorithms of Spectrophotometric Data Correction - Part 11: Use of Splines for 
Approximation of Spectra", IEEE Trans. Instrum. & Meas., Vol. 46, No. 3, June 1997, 
pp. 685-689. 

25 The following operations are accomplished during the initial calibration: 
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• identification of the function g„{^) , using an iterative algorithm as described by 
Jansson; 

• estimation of the parameter ar^ of the fimction g^U) on the basis of a real 

spectrum x'^'^A) using an opiimizaiion algorithm; and 
5 • estimation of coefficients of the rational filter using an optimization algorithm. 

The following operations are accomplished during the continual re-calibration: 

• updating coefficients of the rational filter using an optimization algorithm; and 

• estimating coefficients of the adaptive rational filter using an optimization method 
and a shape of spectrometric data. 

1 0 The exemplary results of spectrophotometric data interpretation, obtained by 

means of the described preferred embodiment are shown in Fig. 6a, Fig. 6b and Fig. 7. 

The following errors are defined to evaluate the exactitude of the procedure AI- 

RC: 

• relative root-raean-square error (RRMSE) 

' It ' ' R 

■ normalized maximum error (NME) 



111 • - NL 

CO 

The errors of estimation of positions and heights of peaks, are: 
^^=0.0015, i?i,= 0.0024; ^i"=0.20, ^^=0.23; 

20 Example 2 

An embodiment of AI-RC was designed for a measuring system which is 
represented by the integrated MM microOSA (Optical Spectrum Analyser) as shown in 
Fig. 8. The analyser is shown used for optical channel monitoring in DWDM networks in 
Fig. 9. 
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The following DWDM system parameters have been selected both for calibration 
and for acquisition of test data: 

• number of optical channels : 8 1 ; 

• wavelength range: A„^;^=1530.77nm, A^,^^=1562.68«w; 

5 • channel spacing : AX = 0.4 nm (50 GHz); 

© laser sources v/ith a known FV/HM = 0.2 rim, where FV.^HM denotes Full Width at 
Half Maximum 

The test data were acquired for a telecommunication bandwidth; its real spectrum 
x{Z) is shown in Fig. 10a. The parameters of this spectrum are the following: 
10 «» the vector of the positions of peaks (carrier frequencies): 
/k= 1530.28 + (k-l) AX + d/k (k=l,..., 81), 
where d/k are realizations of identical independent random variables uniformly distributed 
in the interval [-0.05 nm;0.05 nm], Le. dlj^ oc U (-0.05, 0.05); 

• the vector a of the heights of peaks, that are assumed to be realizations of identical 

15 independent random variables uniformly distribuLi'J in the inten^al [0.01 ; 1], ic, 
fljfcOcU (0.01,L0). 

The adaptation of the MM ^OSA for DWDM applications means that the model 
of the data takes on the form: 

{?«}=Ww**wL.^)+feJ=Ww*^-w*^('^j»»a)L»^)+k ''=1 256 

20 where gj^{X) is the optical response of the spectrometric transducer, 

81 

5(>l;l,a)= ^aj^S{X-lj^) is the so-called idealized spectrum of input light, and [rj^] is a 

sequence modeling the random errors in the data and the amplified spontaneous emission 
(ASE) noise caxised by the laser sources. This is an adequate model for the DWDM 
measurements if the function is an approximation of the shape of the laser light 

25 emission spectrum. 

The a priori information about emission light of a laser is used for reconstruction 
of the idealized light spectrum s{A;l,a). 
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The continual re-calibration sub-procedure allows the measurement of the 
spectrum of laser emission in Is time intervals to continually update parameters p^ of 
reconstruction sub-procedure. 

The chosen operator of projection, mapping the idealized spectrum ^(yl;l,a) into 
5 the space of the data: 

ly„}=G [.(A;I,a);ps ] 
has been defined by the following operations: 

10 y.=yM for /i=i,...,iv 

The function g^{X) has been assumed to have the fonn of the Gauss function: 

1 ( 
gxyi^) = -fr= exp 



Consequently, the vector of the parameters pQ of the operator G contains the discrete 
values of g^(A) , shown in Fig. 1 1 and the parameter cr^ , equal to the spectral 
1 5 bandwidth of the MM |iOS A. 

The chosen operator of reconstruction, transforming the data [y„] into an 
estimate s{X) of j(il;l,a), : 

is defined in the following way: 
20 • a discrete estimate { jc„) of x(X) is found by means of the rational filter applied to 

the data [y^] ; and 

• an estimate s(A) of 5(A;l,a) is computed by means of the spline-based Kalman filter 
applied to {jc„} . 
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The vector = [pb j ,2 • *® parameters of the operator R contains the 
coefficients of the rational filter, as well as the discrete values of the function and 
a regularization parameters required by the spline-based Kalman filter as described in in 

Ben Slima M., Morawski R. Z., Barwicz A.: "Kalman-filier-based Algorithms of 
5 Spectrophctometric Data Correction - Part II: Use of Splines for Approxiir^ution of 
Spectra", IEEE Trans. Jnstinon. cSi Meas., Vol. 46, No. 3, June 1997, pp. 6S5-6S9. 

The following operations are accomplished during the initial calibration: 
• estimation of the coefficients of the rational filter using an optimization algorithm; 
and, 

10 « estimation of the regularization parameter of the spline-based Kalman filter, using an 
optimization algorithm. 

The update of coefficients of the rational filter is achieved during the continual re- 
calibration on the basis of the measured light emission spectrum of a laser. 

Results of adaptive reconstruction and interpretation of teleconinriunicoi-ion data 
15 obtained according to the present exemplar}' embodiment are showTi in Fig. 12. 

The following errors are defined to evaluate the exactitude of the present 
exemplary method: 

■ relative root-mean-square error (RRMSE) 

20 ■ maximum absolute error (MAE) 

A,=maxja„-a„,|}, A,=maxjl„-I^|} 
The errors of estimation of positions and heights of peaks, are: 

£,=210"^%, A,=0.05nm; f^=1.3%, A,=0.5dBm; 
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4. Applicability and scope of the invention 

The proposed procedure AI-RC may be applied in various light spectrum 
measuring instruments and systems, especially in miniaturized integrated spectrometers 

for in situ applications such as environmental probes, optical spectrum analyzers and 
5 ODticn! r^^fc^nnce monitors in telecommmicatinn etc. The motiv?.t!nn for '*t^ 
application in a given measurement situation is founded on expected gai:;s in 
performance and resolution such as the following examples: 

• increased accuracy of spectrometric analyses accomplished by a given spectrometric 
system, resulting from correction of instrumental imperfections and reduced uncertainty 

10 of estimation of the parameters of measured spectra; 

• possibility of automatic calibration of especially miniature integrated spectrometric 
devices during manufacturing process, taking into account their fabrication discrepancies; 

• miniaturization and integration of light-spectrum-measurement related tools, resulting 
from the software compensation of tfieir hardware imperfections and software 

1 5 replacement of some their functions; and 

• reduced costs of spectrometric analyses, given the required accuracy, by replacing a 
high-resolution spectrometric devices with a functionally equivalent but low-resolution 
instrument. 

As is evident to those of skill in the art, a known drawback to miniaturisation of spectral 
20 devices is that manufacturing at very small scales and light propagation through devices 
at those scales often highlights even small imperfections in the manufacturing processes 
and magnifies errors in manufacture. As such, a method of calibrating each device to 
correct for said imperfections would be extremely advantageous in some situations. 

The exemplary embodiment of the invention, presented in section 3 is not 
25 intended to limit the applicability of AI-RC to absorption spectrophotometry and the 
monitoring of optical telecommunication channels. Neither is it intended to limit the 
variety of algorithms that may be used to embody the operations the procedure is 
composed of. On the contrary, the invention is intended to cover alternatives. 
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modifications and. Some practical options for methods embodying the operations of AI- 
RC are briefly characterized. 

4.1. Optional forms of the operators of projection an d reconstruction 

The following mathematical models of the spectrometric data may be used for 
5 defining the operator G : 

a) the stationary linear model: 

b) the non-stationary linear model: 
10 c) the non-linear model, e.g.: 

where giX) and are the apparatus functions of the spectrometric apparatus; 

and Fy are non-linear fimctions. 

15 The corresponding operators G may have the following forms: 

a) the operator corresponding to the stationary linear model: 

V 

where ;,e.. = ^^^)AA; 

b) the operator corresponding to the non-stationary linear model: 

V 

c) the operator corresponding to the exemplary non-linear models: 
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where p^^ .^ - gj^.l^, ^— J 



4.2. Optional methods of signal reconstruction and interpretation 

5 The following methods of signal reconstruction - deconvolution or generalized 

deconvolution - may be used for defming the operator R : 

a) the original domain, numerical differentiation-based method defmed in 
Morawski R. Z., Sokol owski P.: "Application of Numerical Differentiation for 
Measurand Reconstruction", Proc. 7th IMEK0-TC4 Int. Symp, Modern Electrical & 

10 Magnetic Measurements (Prague, Sept. 13-14, 1995), pp. 230-234; 

b) the iterative methods of Jansson and Gold; 

c) the spectram-domain, Tikhonov-regularization-based method; 

d) the cepstrum-domain, Tikhonov-regiilarization-based method defined in 
Kalinowska A., Morawski R. Z., Lubianka T.: "Incorpora- tion of the Positivity 

15 Constraint into a Cepstral Method of Measurand Reconstruction", Proa, Xlll-th IMEKO 
World Congress. (Torino, Italy, Sept, 5-9, 1994). pp. 429-434; 

e) the original-domain, Tikhonov-regularization-based method with the positivity 
constraint imposed on the solution; 

f) the Kalman-filter-based method with with the positivity constraint imposed on the 
20 solution; 

g) the Kahnan-filter-based method with spline-approximation of the solution; 

h) the adjoint-operator method; 

i) the entropy-based variational method; 
j) the Volterra-series-based methods; 

25 k) the rational-filter-based method; and, 

1) the adaptive-rational-filter-based method. 
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Moreover, many other methods developed in the domain of chemometrics as 
described, for example, in Brown S. D., Bear Jr. R. S., Blank T. B.: "Chemometrics", 
Anal. Chem., Vol. 64, No. 12, 1992, pp.22R-49R, and Brown S. D., Sum S. T., Despagne 
F.: "Chemometrics", ^«a/. Chem. Vol. 68, No. 12, 1995, pp. 21R-61R; 

5 telecommunications as described, for example, in Abreu E., Mitra S. K., Marchesani R.: 
"Non-minimum Phase Channel Equalization Usiiig Non-causal Filters", IEEE Trans. 
Signal Processing, Vol. 45, No. 1, Jan. 1997, pp. 1-13 or seismology as described, for 
example, in Berkhout A. J.: Seismic Migration, Elsevier 1985 and image processing as 
described, for example, in numerous publications such as Gonsalves R. A., 

1 0 Nisenson P.:"HST Image Processing: An Overview of Algorithms for Image 
Restoration", Proc. SPIE, Vol. 1567, 1991, pp. 294-307 and Zerwakis M. E., 
Kwon T. M., "On the Application of Robust Functionals in Regularized Image 
Restoration",i'roc/£EE//7r. Conf. Acoustics, Speech & Signal Process. -ICASSP'93 
(Minneapolis MN. USA. April 27-30. 1993)yo\. 5, pp. 289-292 may be applied for this 

15 piupose. 

Tne following methods ars useful for determining the regularir.aiion parameters 
of the operator R : 

a) tlie discrepancy principle with a pre-estiiiiate of the variance of mensurenicnt errors in 
the data as described in Tikhonov A. N., Goncharsky A. V., Stepanov V. V., 

20 Yagola A. G.: Numerical Methods for the Solution of Ill-Posed Problems, Kluwer 1995; 

b) the method of the L-curve as described in Hansen P. C, O'Leary D. P.: "The Use of 
the L-curve in the Regularization of the Discrete Ill-posed Problems", SIAMJ. Sci. 
Comput..Yo\. 14, No. 6, 1993, pp. 1487-1503; 

c) the method of additional set of calibration data as described in Szczecin ski L., 
25 Morawski R. Z., Barwicz A. : "Numerical Correction of Spectrometric Data Using a 

Bilinear Operator of Measurand Reconstruction", Proc. lEEEInstrum. & Meas. Technol. 
Conf - IMTC'95 (Boston. MA. April 24-26. 1995), pp. 488-491. 

4.3. Methods of calibration 

It is often desirable to assume that the isolated peak v,(A,/) has the following 

30 forms: 
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a) the Dirac distribution S{X) for all values of / ; 

b) a triangle whose width is constant or varying versus / ; 

c) a rectangle whose width is constant or varying versus / ; 

d) a Gauss function whose width is constant or varying versus / ; and 
5 e) a Lorenz function whose width is constant or varying versus / . 

The following methods are usefjl in estim.ation of the nppsrnMf; f^,inction g{.?.) : 

a) smoothing approximation applied directly to the data {y^"^ } if the isolated peak 
v,(A,/) is assumed to have the form the Dirac distribution S{X) ; 

b) deconvolutionofthedata {y;"^) with respect to j(A;r^^a'^') ; and 
10 c) subsequent use of deconvolution and smoothing approximation. 

The following methods may be used for determining other parameters of the 
operator R : 

a) a direct transformation of the parameters of the operator G ; 

b) the minimization of any norm of the solution jp^ |j under constraints imposed on 

15 another norm of the disrepancy ||^(A;l'"^a^^')-R [jyr');PR J > ^^'^'^ 

c) the minimization ofany norm of the disrepancy Isl/^r""',:!^^')-!^ [{j^r'jiPR ]|| ^^^^^dcr 
constraints imposed on another norm of the solution ||pr || . 



4.4 Methods for estimation of the parameters of peaks 



The following methods are useful for estimation of the magnitudes a of the peaks, 
20 given the estimates 1 of their positions I : 

a = arg.inf|||{y„}-S [iMl^Ps ]\\^^^} 



or 



i = arg.inf|||5lA)-5(A;i,a)^^ a €a| 
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with A - being a set of feasible solutions; options: q=2 and A c R ; q=oo and A 
c R*;q=2andAc R*;q=ooandAe R*. Some examples of algorithmic solutions 
are given in Morawski R. Z., Mi? kina A., Barwicz A.: "Combined Use of Tikhonov 

Deconvolution and Curve Fitting for Spectrogram Interpretation", Instrum, Science & 
5 Techiohg)^, Vol, 24, No. 3, August 1996, pp. 155-167 and Mr-^y^\ R. Z , Mie kina A., 
Brrv^cz A.: "The Use of Deconvolution nnd Iterative Optimi7?tion S^pectrogrom 
Interpretation", IEEE Trans. Instrum. & Meas., Vol 47, No. 5, Oct. 1997. A particularly 
effective solution of the above optimization problem may be based on a non-stationary 
Kalman filter or an adaptive LMS algorithm is given in Ben Slima M., Morawski R. Z., 
10 Barwicz A.: "Kalman-filter-based Algorithms of Spectrophotometric Data Correction - 
Part II: Use of Splines for Approximation of Spectra", IEEE Trans, Instrum, & Meas,, 
Vol. 46, No. 3, June 1997, pp. 685-689. 

Optionally, the methods for estimation of the magnitudes a are used for iterative 
correction of the estimates of the magnitudes a and positions 1 of the 
1 5 peaks irsterchangeably with the following methods: 



with L being a set of feasible solutions; options: q=2 and L e R * ; q=oo and L c R * ; 
20 q=2andLc Rj;q=ooandLc R*. 

4.5. Mefliods for pre-processing of the spectrometric data 

The following methods are useful for normalization of the data: 
a) the linear or nonlinear transformation of the X -axis, aimed at diminishing the non- 
stationarity effects in the data; 
25 b) the Imear or nonlinear transformation of the ;;-axis, aimed at diminishing the non- 
linearity effects in the data; and 




or 
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c) the linear or nonlinear transformation of the A-axis and j;-axis, aimed at diminishing 
the non-stationarity and non-linearity effects in the data. 

The following methods are useful for smoothing the data: 
a) the linear, FIR-type or IIR-type, filtering; 
i ih- ills- Jim ilUcring; 

c) ihe smoothing approximation by cubic spiiiies; and 

d) the deconvolution with respect to an identity operator. 

Standard methods for baseline correction useful with the present invention are 
described in Brame E. G., Grasselli J., Infrared and Raman Spectroscopy, Marcel Dekker 

10 1976 and in ASTM 1987. (American Society for Testing and Materials): Annual Book of 
ASTMStandards 1987. 

Numerous other embodiments may be envisaged without departing from the 
spirit or scope of the invention, 

15 
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What is claimed is: 

L A method of adaptive inteq^retation of spectrometric data combined with re-caHbration 
for use with a spectrometric apparatus including transducer elements and comprising the 
steps of: 

5 a) initially caiibrating ihc specLrometric apparaais providiiig spcctruiuctiic data; 

b) at intervals automatically re-calibrating the apparatus ab:^enl physical tuning of 
individual transducer elements in isolation by the steps of: 

• providing light from a reference light source to the apparatus, 

• acquiring auxiliary reference data relating to the detected light from the reference light 
10 source, and 

• computing updated results of calibration based on the auxiliary reference data; 

c) sensing a spectrum of Ught {xj with the spectrometric apparatus to provide 
spectrometric data {j^„] representative of the spectrum of light {xj\ 

d) adaptively determining an estimate {s„} of an idealized spectrum {s„} related to the 
15 spectrum of Ught {x^} based on the spectrometric data {v„} and updated results of 

calibration; and, 

e) adaptively estimating positions IxJi^-- and magnitudes 6^i,^2»- of peaks composing 
the spectrum of light {x„ } on the basis of the estimate {f „ } . 

20 1. The method of claim 1 wherein the step (b) of recalibrating involves updating the 

results of a previous calibration absent physical modification and adjustment of the sensor 
elements. 

3. The method of claim 1 wherein the step (a) comprises the steps of: 
25 al) selecting a form of peaks within the idealized spectrum {s„} ; 

a2) choosing the form of the operator R for reconstruction of the idealized spectrum: 
{i^ } = R [{y^}; pj^ ], where p,^ is a vector of the parameters, the form selected in 
dependence upon the selected form of the peaks; 
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a3) choosing rules for adaptation of the parameters to selected features F|, Fj, ... of 

the data {y„ } on the basis of the library of parameters \p^^ \ J = I,...,*/} forming part of 

the calibration data determined during the step of initially calibrating and for updating 
during the step of re-calibrating; 



a8) approximating a relationship between the parameters and the features , , ... 
controlling adaptation. 

4. Tlie method of claim 3 wherein the step (a 6) coiriprises ihc step of: 

15 acquiring a priori information relating to a pre-estimate of the variance of errors in the 
calibration data. 

5. The method of claim 4 wherein the step (a6) comprises the steps of: 
performing elimination of outliers, 

20 performing baseline subtraction, 
performing smoothing, and 
normalizing the data. 

6. The method of claim 3 wherein the step (al) comprises the step of: 

25 selectmg a method of idealizing a spectrum indicative of a form of the peaks. 

7. The method of claim 1 wherein the step (b) comprises the steps: 



:i '] r}Arr-:\r,z n reference liglit nignai for Lhe slcf 





10 {xfO))and{yfO)}fory=l,...,J; 
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bl) choosing a reference light signals for re-calibration: {x^'"'C/)} and related 
for 7=1,... ,7; 

hi) acquiring data } representative of {x;;""' (y)) for J=l,...,J; 
b3) pre-processing the acquired data 

5 b4) updating of the vectors of parameters p^^' on the basis of triplets: [J )]> 
{cr'U)] and tr'O)} for J=U.,J ; 

b5) updating the approximation of the relationship between the parameters P|^ and the 
features Fj , controlling the adaptation. 

10 8. The method of claim 7 wherein the step (b3) comprises the step of: 

acquiring a priori information relating to a pre-estimate of the variance of errors in the 
acquired data. 

9. The method of claim 8 wherein the step (b3) comprises the steps of: 
1 5 perfonning elimination of outliers, 

performing baseline subtraction, 
performing smoothing, and 
normalizing the data. 

20 10. The method of claim 1 wherein the step (d) comprises the steps of: 
dl) pre-processing of [y„ } ; 

d2) computing an estimate of the idealized spectrum {s„} related to the spectrum imder 
study [x„} on the basis of the data {y„] , by means of the operator R with the 
parameters p^^ adapted, at each consecutive step of data processing, to the features 
25 F, , Fj , ... - on the basis of the library of parameters, determined during initial calibration 
and updated during continual re-calibration. 

11. The method of claim 10 wherein the step (dl) comprises the step of: 
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acquiring a priori information relating to accurate estimation of peak parameters. 

12. The method of claim 1 1 wherein the step of acquiring a priori information comprises 
the step of: 

acquiring a priori information relating to a pre-estimate of the variance of errors in the 

<pectroTr:etric data. 

13. The method of claim 12 wherein the step (dl) comprises the steps of: 
perforaiing elimination of outliers, 

performing baseline subtraction, 
performing smoothing, and 
normalizing the data. 

14. The method of claim 1 wherein the step (e) comprises the steps of: 
el) pre-processing of the estimate obtained in step (c); 

e2) computing estimates of the positions and the mngnitisdes of the peaks within the first 
spectrum {:c„}. 

15. The method of claim 14 wherein the estimated positions and magnitudes relate to 
positions and magnitudes within an idealised first spectrum based on the acquired 
spectrum, the updated calibration data, and operators derived therefrom. 

16. The method of claim 15 wherein the step of (el) comprises the step of: 

acquiring a priori information relating to a number of peaks within the spectrometric data 
and an estimate of their widths. 

17. The method of claim 16 wherein the step (el) comprises the steps of: 
performing elimination of outliers, 

performing baseline subtraction, 
performing smoothing, and 
normalizing the data. 



SUBSTITUTE SHEET (RULE 26) 



wo 00/71980 



PCT/CAOO/00583 



32 



18, A spectrometric apparatus comprising: 

transducer elements for sensing a spectrum of light {x„} with the spectrometric apparatus 
to provide spectrometric data {y„} representative of the spectrum of light {x„]\ 
5 means for initially calibrating the spectrometric apparatus in response to a provided 
rcierence light to provide initial results of calibration; 

means for automatically re-calibrating the apparatus in response to a provided reference 
light absent physical tuning of the transducer elements; and, 

a processor for determining an estimate {s„} of an idealized spectrum {s„} related to the 
10 first spectrum of light {jc„} based on the spectrometric data (yj and updated results of 
calibration and for estimating positions IxJi^- and magnitudes a^,a2,^.. of peaks within 
the idealized spectrum {s^] based on the estimate . 

19, A spectrometric apparatus comprising: 

15 transducer elements for sensing a spectrum of light [x^ ] with the spectrometric npparat^js 
to provide spectrometric data {j^^} representative of the spectrum of light {x^] \ and, 
means for initially calibrating the spectrometric apparatus in response to a provided 
reference light to provide initial results of calibration; 
characterised by: 

20 means for automatically re-calibrating the apparatus in response to a provided reference 
light absent physical tuning of the transducer elements 

wherein during calibration of the spectrometric apparatus determination of parameters for 
use in reconstructing an idealised spectrum from data relating to light provided to the 
apparatus is performed. 

25 
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